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Abstract 

Multistable processes are tangent at each point to a stable process, but where the 
^'. index of stability and the index of localisability varies along the path. In this 

n , ! work, we give two estimators of the stability and the localisability functions, and we 

_; ' prove the consistency of those two estimators. We illustrate these convergences with 

-y ■ two classical examples, the Levy multistable process and the Linear Multifractional 

r^ ! Multistable Motion. 
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1 Introduction 
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^ . Multifractional multistable processes have been recently introduced as models for phenom- 

ena where the regularity and the intensity of jumps are non constant, and particularly 
when the increments of the observed trajectories are not stationary. In Figure [H, we dis- 
play a path of a financial data from federal funds, which jumps are more or less marked. 
The multistable processes then extend the stable models in order to take account this 
additional variability (see Figure [2] for an example of a realization of such a process, com- 
puted with the simulation method explained in [4]). We describe then some events with 
a low intensity of jumps at some times, which may be very erratic at other times. We 
provide an other example of application in Figure [7] of Section 14.31 where we consider a 
path coming from ECG data. 

Multistable processes are stochastic processes which are locally stable, but where the 
index of stability varies with "time" , and therefore is a function. They were constructed 
in [H El El [8] using respectively moving averages, sums over Poisson processes, multistable 
measures, and the Ferguson-Klass-LePage series representation, this last definition being 
the representation used here after. These processes are, under general assumptions, local- 
isable, that is they are locally self-similar, with an index of self-similarity which is also a 
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Figure 1: Financial data where the increments do not appear to be stationary : the 
intensity of jumps is varying over time. 




Figure 2: Reahzation of a simulated multistable process. The sample size is n = 20000. 



function. The aim of this paper is to introduce an estimator for each function, the index 
of stability and the self-similarity function. 

Let us recall the definition of a localisable process [21 [3]: F = {Y(t) : t G R} is said to 

be localisable at u if there exists an h{u) G R and a non-trivial limiting process Y^ such 

that 

Y{u + rt) - Y{u) 



lim 



^h{u) 



y:it), 



:i.i) 



where the convergence is in finite dimensional distributions. When the limit exits, Y^ = 
{Y^(t) : t G R} is termed the local form or tangent process of Y at u, and when the 
convergence is in distribution, the process is called strongly h-localisable. 

Ferguson-Klass-LePage series representation 

We define now the multistable processes using the Ferguson-Klass-LePage series repre- 
sentation, that are defined as "diagonals" of random fields that we described below. In 
the sequel, {E,S,m) will be a measure space, and U an open interval of R. We will 
assume that m is a finite measure as well as a cr-finite measure. Let a be a C^ function 
defined on U and ranging in [c,d\ C (0,2). Let f(t,u,.) be a family of functions such 
that, for all {t,u) G t/^, f{t,u, .) G J^a{u){E,S,m). We define also r : E -^ M+ such that 
m{dx) = ■;4^m{dx) is a probability measure. (rj)j>i will be a sequence of arrival times of 
a Poisson process with unit arrival time and (7j)i>i a sequence of i.i.d. random variables 
with distribution P(7j = 1) = P(7j = —1) = 1/2. Let (V^)j>i a sequence of i.i.d. ran- 
dom variables with distribution in on E and we assume that the three sequences (rj)j>i, 
{Vi)i>i, and (7i)i>i are independent. As in [8], we will consider the following random field: 



X{t,u) = Cl{:^^J2^,T;'/''^''\{V,)'/'^(-^f{t,u,V,), (1.2) 
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where C^ = [J^ x~^ sm(x)dx) 

Note that when the function a is constant, then ( 11. 2p is just the Ferguson - Klass - 
LePage series representation of a stable random variable (see [H [TJ [TOl [HI El] and [T5| 
Theorem 3.10.1] for specific properties of this representation). 

Multistable processes 

Multistable processes are obtained by taking diagonals on X defined in (II. 2p . i.e. 

Y{t)=Xit,t). (1.3) 

Indeed, as shown in Theorems 3.3 and 4.5 of [8J, provided some conditions are satisfied 
both by X and by the function /, Y will be a localisable process whose local form is a 
stable process. 

The aim of this work is to estimate both functions of stability a and localisability h. 
Given one trajectory of a multistable process, we provide an estimator for each function 
and we obtain the convergence in all the spaces L^ for the two estimators. We illustrate 
these convergences with two classical examples, the Levy multistable process and the 
Linear Multifractional Multistable Motion. 

2 Construction of the estimators 



Let y be a multistable process defined in (II. 3p . The estimation of the localisability 
function and the stability function is based on the increments (Vfc,Af) of Y. Define the 
sequence (Yk,N)k&,Nen by 

n.. = y(^) - r(^). 

Let to G M fixed. We introduce an estimator of Hito) with 

HN{to) = — — V loglFfc.ivl 

where (?t.(A^))7vgn is a sequence taking even integer values. We expect the sequence 
{HN{tQ))fyf to converge to H{tQ) thanks to the localisability of the process Y. For the 

integers k and A^ such that -^ is close to Iq, , i.l,,. is asymptotically distributed as 

F/ (1). We have then -^ — = HUq) + - — '-— where (ZkN)kN converge weakly to 

logN logN 

k 
— log |1^'^(1)| when A^ tends to infinity and — tends to to. We regulate the sequence 

[Nto] + ^-l 

(Z^^n) near to using the mean ^^jW ^ Z^^jy and we can expect this sum will be 

k=[Nto]-^ 



bounded in the U spaces to obtain the convergence with a rate j^^- The convergence is 
proved in Theorem 13.11 

Let po > and 7 G (0, 1). With the increments of the process, we define the empirical 
moments Sn{p) by 

- y 



SNip) 

Let 



n{i\) ^-^ 

\ fc=[jvto]-^ / 



^ and fi„(rt ^ '^i!!:;;;" 



^exp(p) = ^TTT ^^^ ^"(P) = /irlVl^M/^ ^P<" 



where Z is a standard symmetric a-stable random variable (written Z ~ ^^(l, 0, 0) as in 
^Lojj, i.e t|^| - pj+oo^_p_i^.^2(„)^„ • 

Consider the set A^ =: argmin ( J |i?exp(p) — -Ra(p)r'^p) • Since the function a — > 
/ l-Rcxp(p) — -Ra(p)r'^P 1 is a continuous function, A^v is a non empty closed set. We 

\P0 J 

define then an estimator of a (to) by 



On (to) = min argmin ( / |-Rexp(p) - Ra{p)Vdp J 
\ Qefo,2i \Jpo J 



-2 \ 1/7^ 

'(/ 

:Ke[0,2] \Jpo 

Under the conditions of Theorem 13.21 Y is if (to)-localisable and ^^^(l) ~ Sa{to){l, 0, 0) so 

\Y 1^ 

-^ ' , . converge weakly to [^^'^(l)^ and with a meaning effect, N^^*°'^Sn{p) tends to 

{E.\Yl^^{l)\^Y''^ in probability, which is the result of Theorem l3.2[ Without more conditions, 

/p^ |i?exp(p) - Ra{p)ydp tends to /^^^ \Ra{to){p) - Raip)\"'dp . Naturally, a(to) is the only 

solution of argmin^g[Q2] / \Ra(to){p) — Ra{p)\'^ dp and this leads to the definition of dAr(to)- 
The convergence is proved in Theorem 13.31 

3 Main results 

The three following theorems apply to a diagonal process Y defined from the field X 
given by (IL2I) . For convenience, the conditions required on X and the function / that 
appears in (11.21) . denoted (CI), . . . , (C14), are gathered in SectionEl Theorem 13. II lead to 
the convergence in the U spaces of the estimator of the localisability function H, while 
the two Theorems 13.21 and 13.31 draw to the convergence of the estimator of the stability 
function a. 

3.1 Approximation of the localisability function 

Theorem 3.1 Let Y a multistahle process. Assume the conditions (CI), (C2), (C3) (or 
(CI), (Cs2), (Cs3) and (Cs4) in the a-finite space case), and that there exists a function 
H such that (C5)-(C14) hold. Assume in addition that lim -7^ = +00. 



N^+00 



n{N) 



Then, for all to E U and all r > 0, 

lim E H^ito) - H{t,] 

Proof 

See Section [5l 

Remark: Under the conditions (CI), (C2), (C3) and (C5) listed in the theorem, 
Theorems 3.3 and 4.5 of [8] imply that Y is if (to)— localisable at to- 

3.2 Approximation of the stability function 

We first give conditions for the convergence in probability of Sn{p) in Theorem 13. 2[ which 
is useful to establish the consistency of the estimator dAr(to). 

Theorem 3.2 Let Y a multistable process. Assume the conditions (CI), (C2), (C3) (or 
(CI), (Cs2), (Cs3) and (Cs4) in the a-finite space case). Assume in addition that: 

• lim n(N) = +oo. 

N->-+oo 

• lim -7^ = +00. 

• The process X(.,to) is H [to) -self- similar with stationary increments and i/(to) < 1. 

• (C*) There exists ei > and jo £ N such that for all j > jo, 

ho^t^{x)hj,t,{x)l\^^m{dx) < (1 - ei)llVtolla(Jo)' 

E 

where hj^u{x) = f{j + 1, u, x) - f{j, u, x). 

• lim fj^\ho^to{x)hj^to{x)\^^m{dx) = 0. 

Then, for all p G \po, a{to)), 

Ar^(*«)5jv(p) ^ {E\Z\Py/P 

Af— >-+oo 

where the convergence is in probability and Z ~ Sa{tQ){i, 0, 0). 

Proof 

See Section O 

Theorem 3.3 Let Y a multistable process. Assume the conditions of Theorem \3.2\ then, 
for all to E U and r > 0, 

lim E |djv(to) — a(to)r = 0. 

Af— >-+oo 

Proof 

See Section [5l 



4 Examples and simulations 

In this section, we consider the "multistable versions" of some classical processes: the a- 
stable Levy motion and the Linear Fractional Stable Motion. We provide then an example 
of application with ECG data. 

We first recall some definitions. In the sequel, M will denote a symmetric a-stable 
(0 < a < 2) random measure on R with control measure Lebesgue measure L. We will 
write 

L„(t) := /" M{dz) 
Jo 

for a-stable Levy motion, and we will use the Ferguson-Klass-LePage representation, 

oo 

vt G (0, 1), L„(t) = cy- j]7.rr'/"i[o,*](^.). 

The following process is called linear fractional a-stable motion: 

/oo 
fa,H{b^,b~,t,x)M{dx) 
-oo 

where t eR, H E (0,1), b+, b' E R, and 

/,,^(6+, b-,t, x) = 6+ ((t - x)f -^/° - (-X) J-^/°) 

+ b- ((t - x)^-^/" - (-a:)^-^/" 

When fe"*" = 6^ = 1, this process is called well-balanced linear fractional a-stable 
motion and denoted La^n- 

The localisability of Levy motion and linear fractional a-stable motion simply stems 
from the fact that they are 1/a-self-similar with stationary increments [3]. 

We now apply our results to the multistable versions of these processes, that were 
defined in [H [5] . 

4.1 Symmetric multistable Levy motion 

Let a : [0, 1] — )■ [c,d] C (1,2) be continuously differentiable. Define 

oo 

Xit,u) = C^iCJ") 5^7.r-^/"(")l[o,](V.) (4.4) 

and the symmetric multistable Levy motion 

Y{t)=X{t,t). 
Proposition 4.4 If lim n(N) = +cxd and lim -t^ = +oo, then for all r > 0, 



lim E 



HNito] 



a (to 



and lim E |djv(to) — a(to)r = 0- 

N-^+oo 



Proof 

We know from [9] that all the conditions (C1)-(C14) are satisfied. We deduce from 

Theorem 13.11 that lim E HN(to) ttt = 0. Since the process X(.,to) is a Levy 

motion a(to)-stable, X(.,to) is ^^Tpr-self-similar with stationary increments [15]. We then 
prove that the condition (C*) is satisfied. 



so for j > 1, 

a(to) 

ho^t^{x)hj^t^Xx)\ 2 dx = 0. 

R 

We conclude with Theorem 13. 31 b 

We display on Figure [3] some examples of estimations for various functions a, the 
function H satisfying the relation H{t) = -^. The trajectories have been simulated 
using the field (14. 4p . For each u G (0, 1), X{.,u) is a Q;(M)-stable Levy Motion. It is then 
an a(M)-stable process with independent increments. We have generated these increments 
using the RSTAB program available in [16] or in [15], and then taken the diagonal X{t, t). 

Each function is pretty well-evaluated. We are able to recreate with the estimators 
the shape of the functions. However, we notice a significant bias on Figure El in the 
estimation oi H. It seems to decrease when H is getting values close to 1. We observe 
this phenomenon with most trajectories, while the estimator a seems to be unbiased. We 
have displayed the product aH in order to show the link between the estimators. We 
actually find again the asymtpotic relationship H{t) = ^4^. 

We observe on Figure H] an evolution of the variance in the estimation of a. It seems 
to increase when the function a is decreasing, and we conjecture that the variance at 
the point to depends on the value a (to) in this way. In fact, the increments Yk^N are 
asymptotically distributed as an a(to)-stable variable, so we expect that Sn and i?exp 
have a variance increasing when a is decreasing. 

We have increased the resolution on Figure taking more points for the discretization. 
The distance observed on Figure |Hb for a near 1 is then corrected. 




0.2 3 4 0.5 



a{t) = 1.98 - 0.96t 




0.2 3 4 0.5 



a{t) = 1.98 - 



0.96 



1 + exp(20 - 40i) 





0.2 3 4 0.5 



Hit): 



1 



1.98 - 0.96t 

Eslimilion oT h 




0.2 3 4 0.5 



Hit) 



1 + exp(20 - 40t) 
1.02 + 1.98 cxp(20-40i) 



0.2 3 4 0.5 





0.2 3 4 0.5 




0.2 3 4 0.5 




0.2 3 4 0.5 



ait) = 1.5- 0.48 sin{27rt) 



Hit) 



1.5- 0.48 sin(27ri) 



Figure 3: Trajectories on (0, 1) with A^ = 20000 points, n{N) = 2042 points for the 
estimator a, and n{N) = 500 for H. a and a are represented in the first column, H and 
H in the second column, and in the last column, we have drawn the product aH. 




0.2 3 4 0.5 





0.2 3 4 0.5 



a) b) c) 

Figure 4: Trajectory of a Levy process with a{t) = 1.5 + 0.48sin(27rt) in figure a), and 
the corresponding estimation of a in figure b) with n{N) = 2042. The figure c) represents 
various estimations of a for the same function a{t) = 1.5 + 0.48sin(27rt), with different 
trajectories. 



8 



i\ff 


k/1 


r^ 


\V 



0.2 3 4 0.5 




0.2 3 4 0.5 



d) 



e) 



Figure 5: Trajectory with A^ = 200000 in figure d), and the estimation with n{N) = 3546 
in figure e). 



4.2 Linear multistable mult ifr actional motion 

Let a : R — )■ [c, (i] C (0, 2) and iJ : R — )■ (0, 1) be continuously differentiable. Define 

2-2 °° 

(4.5) 
and the linear multistable multifractional motion 

Y{t)=X{t,t). 



Proposition 4.5 Assume that H — - is a non-negative function, lim n(N) = +oo and 



lim ^^7^ = +CXD. r/ien for all r > 0, 



Af-s>+oo 



lim E 

N-^ + QO 



-f^Af(^o) ~ -f^li^o) = and lim E \aN(to) — Ci(to)\^ = 0. 



Af- 



Proof 

We know from [9j that all the conditions (C1)-(C14) are satisfied. We deduce from The- 

orem l3.1l that lim E i^Ar(to) — H^to) = 0. Since the process X(.,to) is a {H(to),a(to)) 
linear fractional stable motion, X(.,to) is i/(to)-self-similar with stationary increments 
[T5] . We write /ij,to(a;) = (^(j— x) with (yf(-u) = |l+-u| ° "Co) —\u\ "^"^Co) so with Propo- 



sition 2.2 of [T3]. the condition (C*) is satisfied. Let us show that lim /r |^o,to(^)^j,to(^)l ^ ^^ 



and conclude with Theorem | 

Let e > 0. Let Cq > such that / \hQ^toi^)\"^^°^ dx < |. Since Vx G [— co,co] 

|x|>co 

lim |/io,io(a^)^i,to(^)l^^ = ^^d {hj^to{^))j is uniformly borned on [— Co,Co], 



j- 



lim / |/io,io(a;)/ii,to(a;)|^~<ia; = 

^^+°°J\x\<ca 



'\x\<co 

We show on Figure [6] some paths of Lmmm, with the two corresponding estimations 
of a and H. To simulate the trajectories, we have used the field (14. 5p . All the increments 
of X{.,u) are {H{u),a{u))-\meaY fractional stable motions, generated using the LFSN 
program of [T6j. After we have taken the diagonal process X{t,t). 

These estimates are overall further than the estimates in the case of the Levy process, 
because of greater correlations between the increments of the process. However, the 
estimation of H does not seem to be disturbed by those correlations. The shape of the 
function H is kept. For a, we notice some disruptions when the function is close to 1. We 
finally show an example where the estimation of a is not good enough in the last line of 
Figure [61 The trajectory. Figure El a), seems to have a big jump, which leads to decrease 
the estimator a, represented on Figure |6lb), while the jump is taken account in the n{N) 
points. The estimation of H, represented on Figure Oc), does not seem to be affected by 
this phenomenon. 
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ait) = 1.41 + 0.57i 



H{t) = 0.725 + 0.175 sm(27rt) 






a{t) = 1.695 + 0.235 sm(27rt) 



H{t) = 0.725 - 0.175 sm(27rt) 








a{t) = 1.695 + 0.235 sin(27rt) 



Hit) =0.59 + 0.31* 



Estimalionofdplini 






a{t) = 1.41 + 



0.47 



1 + cxp(20 - 40i) 





Hit) =0.9-0.35* 



„= 


rAp'MLi 


Q85 
0.75 


1 / . 







a) 



b) 



c) 



Figure 6: Trajectories with N = 20000 in the first column, the estimations of a with 
n{N) = 3000 points in the second column, and in the last one, the estimations of H with 
uIn) = 500 points. 
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4.3 Application to ECG data 

We consider an example of trajectory with a varying index of stability and a varying index 
of localisability. The dataset conies from [12] . 

We denote Z the process corresponding to an ECG records. Its length is A^ = 1000000 
points. We consider then the process Y defined by 



J-'O) = E ^w 



fc=i / 



The realization of process Y associated to EGG series is represented in Figure [71 The 
increments of this process can not be regarded as stationary. We see in this example that 
the smoothness, as the intensity of significant jumps, is actually varying with time. 



^Cl:C3 dsLis 



-O.^l - 




Figure 7: Trajectory of the process Y associated to EGG series with A^ = 1000000 points. 

We have done an estimation of the localisability function H for this process Y . Figure 
[8] represents an estimation of H as function of t. The estimate of H is calculated by taking 
n{N) = 25000 points. 



t i nri^t icTi 1-1 c^-f I 



□ T^e - 




Figure 8: Estimation of H calculated on the process represented in Figure [71 

We notice a correlation between the noisy areas of the trajectory and the times when 
the exponent H is small, and also a greatest exponent when the trajectory seems to be 
smoother. For the estimation of the function a, we have taken ?t,(A^) = 25000 too. The 
result is presented in Figure M We observe also here a link between the noise and the 
function a. When the intensity of the significant jumps of the trajectory is high, the 
stability function is close to 2. A lower stability index matches to a period with a lower 
intensity of significant jumps. 
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o.-i o.s 0.3 o.^ o.s o.e o.^:^ o.s o.s 



Figure 9: Estimation of a calculated on the process of Figure [71 



5 Proofs 



Proof of Theorem 13.11 

Note that it is sufficient to prove the result of Theorem 13.11 for r > 1 since the 
convergence in L^ implies the convergence in L'^ for all q < p. Let r > 1. Let H satisfying 
the condition (C5). We write 



H^ito) - H{to) 



[A'tol + ^-l 



n(A^)logA^ ^-^ 



log 



Yi 



k,N 



l)Hito) 



N 



n(N) hgN liNtp] njN) 



[Ntp] n(N) 
N ~'~ 2N 



log 



Y{ 


lNt] + l. 
N > 


-y{ 


[NtU 
N ) 




(^)H(.o) 





dt. 



Let 6N{dt) 



N 



l|- [iVto] n{N) ^. [Ntp] n(N) -,dt aud f N {t) — lOg 



"V^V i jv 2N -''^ N '*' 2N } 

Since Jg 6iy{dt) = 1, we obtain 



n^^^)-y(^) 



(i)^w 



HN{to)-Hito) 



, ,, , fN{t)6N{dt)+ {H{t)~H{to))6N{dt). 
log iV Jo io 



Then, there exists a constant Kj. E R depending on r such that 



l-f^Afli^o) — H{tQ 



E(\J,'h{t)6N{dt)\^ 
< Kr^ — '- + K, 



\ogN\ 



{H{t)-H{t,))5N{dt) 



H is continuously differentiable and lim 



N 



N~^+oo "(^) 



'0 

+00 so 



lim I {H{t) - H{to)) 6N{dt) = 0. 



To conclude, it is sufficient to show that there exists a constant K E H depending on 
to and r such that for all A^ G N, E f| J^ fNit)5N{dt)\^] < K. Let U an open interval 
satisfying all the conditions (C1)-(C14), and to G U. We can fix A^o ^ N and V d U an 
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open interval depending on to such that for all N > Nq and all t G V", * ^ G U, -LJ- G U 
and Jq /Af(t)5jv((it) = Jy fN(t)SN{dt). With the Jensen inequality, 



/jv(t)5^(rft)r < / E\Mt)r5^{dt). 



We consider E|/jv(t)|'' = / P (|/jv(t)r > a;) rfa;. 





+00 



Ei/^(t)r 



P (/^(t) > x^/") dx+ P (/jv(t) < -x^/") da; 



nM±l)_w[^^]' 



A^ 



A^ 



,1/r 



> 



Ar-f^(t) 



dx 



+ CO 









wM±l)_w[^^]' 



AT 



N 



< _„.,^ I dX. 



N^(t) j 



y,m + l [Nt]. 



N 



N 



,1/r 



> 



N"(t) 



< P 



y,M + i)_KW 



+ P 



A^ 



r(t)-ri 



a/r 



> 






2A^^W 



> 



2A^^(*) 



so 



with 



and 



Ei/^(t)r<i]v(t)+i^(t)+i^(t) 



+00 

l]vW = / P 



+00 

i^W = / p 





„M±i)_y(,) 



Fi 



iV 
■ A^ 



> „,., I ax. 



r(t) 



> 



2Ar^wJ 



dx 



l%{t) 



Y\ 



[Nt] + 1^ ^^iNt] 



N 



Y( 



N 



,1/r 



< 



Ar^(t) 



dx. 
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We consider first Ijv(^)- 



+00 



+00 





A-(M±l,M±i)^A-(M±i.() 



Xi 



N ' A^ 
[Nt] + 1 



A^ 



,1/r 



> 



iN^^^^i 



dx 



N 



,t)-X{t,t) 



> 



4iV^(t)J 



dx. 



With the conditions (CI), (C2) and (C3) (or (CI), (Cs2), (Cs3) and (Cs4) in the 
cr-finite space case) we can apply Proposition 4.9 or 4.10 of [9] : there exists Ku > such 
that for all (m, v) G f/^ and a; > 0, 



P{\X{v,v)-X{v,u)\>x)<Kc 



V — u\ 



X" 



:i + 1 loe 



\v — u\ 



X 



+ 



\v — u\ 



x"- 



:i + 1 log 



so there exists Kjj > such that for all N > Nq and all t E V, 



xfM±l,M±i)_xfM±l,,) 



+Ku 



N ' N 
(log NY 

^ I JYc(l-H(t))gCa'l/r 

(logiV)'^ 

]^d{l-H{t))^dx^/^ 



N 






+ Kr 



X 



c/r 



+ K, 



u 



JYc(l-_ff(t))gCxi/'- 

^d/r 
]\Jd{l-H{t))gdxi/^ 



\V — U\ 
X 

(5.6) 



Since H^ < 1, we conclude that 



'+ ^ J-5 



lim / sup P 





A-(M±i,M±i)_.«M±i,() 



N 



N 



N 



Air 



> 



4iV^(*) 



dx = 0. 



With the same arguments. 



lim / sup P 





.Y,MM)_.v(M() 



A^ ' A^ 



,,1/r 



Let rj < c. The Markov inequality gives 



.V(^,*)-X(M) 



> 



4iV^(*) 



A^ 



AVMVHit) 
< r^E 



> 



4A^-f^(*) 



t/x = 0. 



(5.7) 



^VX 



l/r 



ix(^^^^,t)-x(t,t)r 



and Property 1.2.17 of [T5 

.[ATij + l 



X( 



N 



-,t)-x(t,t)r 



Ca(t),o(r/)'' ( /^l/(^^^^,^,^) - fit,t,xr^'^) m{dx)^ 



With the condition (C9), there exists K > Q such that for all N > Nq and all t eV, 
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A-(li!^.<)-.Y(M) 






Finally, there exists K > such that for all N > Nq and all t G V^, 

I],{t) < K. 

Using the equation fl5.7p and the condition (C9), we obtain that there exists i^ > such 
that for all N > Nq and allt eV, 

lUt) < K. 

Thanks to the conditions (CI), (C6), (C7), (C8), (CIO), (Cll), (C12), (C13) and 
(C14), we conclude for I^(t) using Proposition 4.11 and 4.8 of [9] : there exists K > Q 
such that for all N > Nq and all t & V, 



[iVt]+^ [iVt] 

^ N ^ ^ N ^ 



< 






a/r 



<K 



Ar-f^(t) 






so there exists K > such that for all N > Nq and all t & V, 

I%it)<K , 

Proof of Theorem [3:2] 

Let p G \po,a(tQ)). We define 



ANip) 



j^,Hito) [^*ol + -H)_, 



n{N) 



E 



k=[Nto]-'-^ 



Xi 



k + l k + 1, 

■ N ' N ' 



Xi 



k + l 



to) 



and 



5iv(p) 



CNip) 



JVP-H'(to) 

n(X) 



[A^to] 



n(iV) 



E 



fe=[7Vio]-^ 



^(f-f)-^-(|-*»' 



]\fpH{to) ^ "' 2 



n{N) 



E 



fc=[ivto]-^ 



^{ jy ) ''^O) - ^(]y' '^O; 



We have, for p < 1, 

P (|X^^(*«)5^(p) - E|Zn > x) < P (|X^^(*«)5^(p) - C^(p)| > ^) + P (|E|Zr - C^(p)| > 

< P (|E|Zr - CAr(p)| > I) + P (^iv(p) + i?iv(p) > 



and for p > 1, 

P f |X^(*«)5iv(p) 



:Eizn^i >x 



< P |X^(*«)5^(p)-C^(p)|> 



X 



p ic^(p)-(E|znF|> 



X 



< p ( i(E|zn^ - cUp)\ > I) + p (4(p) + bUp) > I 
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p p 

To prove Theorem 13.2^ it is enough to show that Ai^{p) — )■ 0, B^ip) — > and 



N 



nfAO'-^-r '^'oi "(^)</^ i^*oi I «(N)^dt. Let U 



C^(p) -^ E|Z|^ with Z ~ 5„(t„)(l,0,0). 

p 
We consider first A]\r{p) — )■ 0. Let Sn^cH) 

an open interval satisfying the conditions of the theorem and to G U. We can fix A^q ^ N 
and V G U an open interval depending on tg such that for all N > Nq and all t G y, 
- — '— G U, ^-^ G ?7, Jg SNidt) = Jy6N{dt), and such that the inequality f l5.6p holds. 



AT 



P(Ajv(p)>x) 



y-^ [Nt] + l [Nt] + 1 ^ __ v/ [Nt] + l 



N ' N 



■)-X(^^,to) 



< 



X 



- / E 

X Jv 



(l/iV)^(*o) 
Yf [^*l + 1 [Nt\ + l^ _ ^f[Nt\ + l 



N ' Af 



lV~'''OJ 



{l/N)Hito) 



X 



■ \Nt\ + \ [Aff]+1 ' 
. Af ' iV ' 



Xi 



[Aff]+1 
Af ' 



to) 



(l/iV)H(to) 



bfq{d£) > X 
6N{dt) 
(5 AT ((it) 



Let t G V". 



N ' Af 



1V~'''0J 



(l/Ar)^(to) 



■ [Aff]+1 [Aft] + l ^ 

, Af ' Af ' 



• [Aft]+1 



(l/Ar)^{to) 



> M 



1/P rfw. 



Let u > 0. We know from (J5.6P that there exists K^/ > such that for all t E V, 



X 



' [Nt]+1 [Ni\ + 1- 
■ Af ' Af . 



X 



' lNt]+l . , 
, Af ' 0^ 



(l/Ar)^(to) 
so, with the assumption H(to) < 1, 

^([Af|tl,[AfjM)_X(MLti,to) 



y\ ((logiV)^+|log^r) ((logiV)'^ + |logn|'^) 



lim P 

Af^-+oo 



(l/iV)^(*o) 



> M 



1/p 



0. 



There exists Ku,p > such that 

^ [ N ' Af ) ^ ^ I Af ''^OJ 



(l/Ar)^(*o) 



> u^^" 1 < l«<i + Ku,p 



logMl"^ I logM 



U 



d/p 



+ 



u 



c/p 



L«>1- 



(5.^ 



Since a is a continuous function, we can fix U small enough such that c = inf a(t) > p. 

teu 

We deduce from the dominated convergence theorem that for all t G U, 



lim E 

Af-i>+oo 



X 



■[Nt\j-1 [Aft] + l ^ 

, Af ' Af ' 



■X(I^,to) 



(l/Ar)^(to) 



0. 



With the inequality (15. Sp . 



X 



■ [Nt]+1 [Nt\+1- 

■ Af ' Af . 



X 



■ [Aff]+1 , ^ 
, Af ' 0^ 



(l/N)^{to) 



< 1 



'"~AVJ^^ + ^^1* 



u 



d/p 



u 



c/p 
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and again with the dominated convergence theorem, 

hm P {An{p) > x) = 0. 

p 

The same inequahties holds with Biy{p) so we obtain B^{p) — > 0. We conclude proving 

p 
Cn{p) — > E\Z\P. Let Co > 0. We use the decomposition 



[A^to] + ^-l 



C^{p)-E\Z\ 



n{N) 



E 



k=[Nto]-^ 



-^(l^i'^o) --^(jv'^o) 



(l/Ar)^(*o) 



^ x i!^,to)-x(-^,to) 



(l/N)H(to) 



>C0 



-^z\n\z\> 



CO 



[Nto] 



n(N) 



n{N) 



E 



fc=[iVto]-^ 



-^(lv~'^o) ~^(lv'^o) 



(l/Ar)^{to) 



1 x(fc+i,to)-x (-fc,,to) 



(l/JV)-f^(*0) 



<C0 



Ei^riizi 



<co- 



Let e > and a; > 0. By Markov's inequality, we have 



/ 



Pi 



[JVto] + 4^-l 



- y 

n(N) ^ 



-^(lv~'^o) ^ ^(tv'^o) 



:i/A^)^(*o) 



■'■ ^(n^.to)-^ (l^.to) 



[Afto] 



n(iV) 



< 



xn{N) 



E E 



fc=[iVto]-^ 



-^(l^5^o) - -'^(tv'^o) 



(l/Ar)^(to) 



(i/iv)-f^(*o) 



■y(Ty^.*o)--y(|-.*o) 



(1/]V)-H"('0) 



a; 

>co 4 



>co 



\ 



Since X(.,to) is -ff(to)-self-similar with stationary increments, 

Pi<^E[|X(l,to)ri|x(Mo)l>co] 
and 



Ei^riizi 



<co 



k E E 



n(A^) 



k=[Nto]-^ 



-^(lV~'^o) - -^(iv'^o) 



(l/Ar)H(to) 



^(T^.to)-^(^io2 



(1/]V)-H"('0) 



<Co 



We fix Co large enough such that for all A^ G N, Pi < | and E\Z\p1\z\>c„ < f • Writing 
K{x) = \x\^l\x\<co and AXk^ = X(A; + l,to) ~X{k,to), using Chebyshev's inequality, we 
get 



P{\Cn{p)-E\Z\p\>x) < J + -3-1^ J2 Cov{K{AXk,t,),K{AX,,,)) 



2 x^n{Ny 



k,j=[Nto]~^ 



n(N)-l 



- 9 + ;:^ ^(AO + ^;^ 4^ CoviKiAXo,,),KiAX,,,,)) 



2 a;2 
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Under the condition (C*), we can apply Theorem 2.1 of [13]: there exists a positive 
constant C such that 

Cov{K{AXo,to),K{AXj,t,)) < C\\K\\l [ \ho,to{x)h,,to{x)\'^m{dx). 

Je 

Since the process X(.,to) is -f/^(to)-self-siniilar with stationary increments, the constant 
C does not depend on k, j. We then obtain the existence of a positive constant Cp^co 
depending on p, cq and x such that 

Since hm n{N) = +00 and hm /"„ |/io,io(a^)^i,to(^)l^^"^('^^) = O5 "^^ conclude 
with Cesaro's theorem that there exists Nq & N such that for all N > No, 

-^y^|Vo(x)r(*«M^x) + ^ -£ j^\ho,,{x)h,,,{x)\'^m{dx) <'- 

and 

P(|C;v(p)-E|Z|P| >x) <e ■ 

Proof of Theorem 13:3] 

Since x — )■ x"^ is an increasing function on R+, 

aAr(to) = min argmin / |-Rcxp(p) - Ra{p)pdp . 

\ ae[0,2] 7po J 

2 2 

Let 5'Ar(a) = / |-Rcxp(p) - Ra{p)Vdp and 5f(a) = / |i?a(io)(p) - Ra{p)Vdp. 

Po pa 

g is a. continuous function on (0, 2], with g{Q) > 0, g{2) > 0. The only solution of the 
equation g{a) = is a(to). Moreover, lim | _~^/^ -.i-^ > 0. 

Then, there exists Ka^to) a positive constant depending only on a(to) such that: 

VaG(0,2), \g{a)\>K^^to)\a-a{to)\. (5.9) 

We estimate now \g{QN(to))\. 

\g{aN{to))\ < \giaNito)) - gNiaNito))\ + \9NiaNito))\ 
< \g{aN{to)) - gN{aN{to))\ + gN{a{to)), 

and 

Igi&Nito)) - gNi&Nito))] 

'po 
2 



{\Ra{to)ip) - Ra^{to)ipW - \RcKpip) - RaN(to)iPW) ^P 

< / \Ra{to) (P) - Re^p (P) I "^ dp 
-'PO 

gN{a{to)). 
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From §M>, 



\aN(to) - a(to)\ < — g{aN{to)) 



K. 



a(to) 



< 



K. 



a{to) 



-gN{a{to))- 



Let us show that Vr > 0, hm E |(7jv(a(to))r = 0- Let r > 0. One has, using the 

A''— >-+oo 

inequahty Sj\f{p) < S]\f{q) for p < q, 

a{to) 2 

9N{a{to)) = / \Re^p{p) - Raito){p)Vdp+ / \Rey,p{p)Vdp 



PO 

Q(to) 



a (to) 



< / \Rc.p{p)-Ra{to){p)\''dp+{2-a{to)) 

PO 



Sn{Po 



S'jv(a(to)) 



For the first term, we use Theorem 13.21 : for all p G [po, a(to)), 

N^'^'^^SNip) -^ (E|Z|P)i/P (5.10) 

where Z ~ S'a(jo)(l, 0, 0). It is clear that Wp G [^05^(^0)), 

(iV^(*«)^jv(po),A^^^*''^57v(p)) ^ ((E|Zr)^/P«,(E|Zn^/P), 
and 

^exp(p) = ^^^A/?.(,„)(p). (5.11) 

Note that VA^ G N, Vp G [po,a{to)), |i?exp(p)| < 1 so there exists a positive constant 
K depending on •yr, a (to) and p such that 



E|i?cxp(p) -i?a(ta)(p)| 



7r 



A' 



P (|i?exp(p) - i2a(t„)(p)r'^ > x) C?X. 



Finally, with (15. lip . Vp G [po? «(^o))> E|i?cxp(p) ^-Ra(to)(p)r^ — ^ 0- With the inequality 

E|-Rexp(p) — Ra{to){p)\'^^' < 2C^r where C^r is a positive constant depending on •yr, by the 
dominating convergence theorem, 

r-Q(to) 



lim / E|i?cxp(p) - Ra{tQ){p)r'' dp = 0. 

^+°%o 



7V-S-+C 



To conclude we show that 

SnJPo) 
5jv(a(to)) 



SnJPo) 



5jv(Q!(to)) 

to show /y,^^r\, -^ 0. Let p < a(to)- 



1 



0. Since ViV G N, 



■5'jv(po) 



iiV^(*o)^^(a(to))| 



> x) < P( 



1 



|iV^(*»)^7v(p)| 



5]v(a(io)) 



> x) 



< 1, it is enough 
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So, 



1 

lim sup P I ^^o-/^ ^ ^ , , — TTT 



> x) < limsup P(--— TTT-TTT^-TT 

1 



> X) 



lim Pf — 



> x) 



1 



with (Eini). Since ^Jim^^ P( (Eiz!p)i/p > 3^) = 0' ^e have hmsup P( |jvi^(.o)g^(^(to)) 



>a;) =0 



and 



7V-»('o)5iv(a(io)) 
giv(po) P 



P - -• - ,.m/„^., / ^ P 'E|Zpo)i/P«, we obtain 



0. Using the convergence N^^^'^^S^^po) 



SN{a{to)) 







6 Assumptions 



This section gathers the various conditions required on the considered processes so that 
our results hold. 



(CI) The family of functions v — ?■ f{t,v,x) is differentiable for all (f,t) in f/^ and 
almost all x in E. The derivatives of / with respect to v are denoted by /^. 



(C2) There exists S > - -1 such that : 



sup 
teu JR 



snp{\f{t,w,x)r^-^) 

well 



1+5 



rh{dx) < oo. 



(Cs2) There exists 6 > - — 1 such that 



sup 

teu Jn 



sup(i/(t,^,x)r('")) 

weu 



ll+S 



r(x) m(dx) < oo. 



(C3) There exists S > j; — 1 such that : 



sup 
tec/ Jn 



sup(|/:(t,u;,x)r("')) 

wGU 



l+S 



m{dx) < cxD. 



(Cs3) There exists 6 > - — 1 such that : 



sup 



sup(i/:(t,^,x)r('")) 



ll+S 



r{x) m{dx) < oo. 



(Cs4) There exists 6 > - - I such that : 



sup 

teu Jn 



sup 

well 



|/(t,w,x)log(r(x))|"^"'^ 



1+5 



r(x) m(dx) < oo. 
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(C5) X{t, u) (as a process in t) is localisable at u with exponent H{u) G (-f/^-, Hj^) C 
(0, 1), with local form X^(t,-u), and u h-> -ff(-u) is a C"*^ function . 

(C6) There exists K[/ > such that \/v eU,\/ue U, \/x E R, 

\f{v,u,x)\ < Ku- 

(C7) There exists Ku > such that \/v eU,\/ue U, Vx G R, 

\f^{v,u,x)\ < Ku- 



(C8) There exists Kjj > and a function /J defined on U such that Vf G f/, Vm G f/, 
Vx G R, 

(C9) There exists eo > 0, Ku > and a function iJ defined on U such that Vr < ^o, 

\/teU, 

-^^ j^\f{t + r,t,x)- f{t,t,x)r^^m{dx)<Ku. 

(CIO) There exists p G (c?, 2), p > 1, -ft'c/ > and a function iJ defined on U such 
that Vt; G f/,VM G f/, 



|/(f,-u,x) — f{u,u,x)\^ m{dx) < Ku- 



I |l+p(_ff(n) i^) ,„ 



(Cll) There exists i^^y > such that Vf E U,Wu E U, 



\f{v,u,x)fm{dx) < Ku. 



R 



(C12) 



^et/ 



inf / f{v,v,xYm{dx)>0. 



R 



(C13) There exists a positive function g and a function iJ defined on U such that 
1 



lim sup 



j,l+2{H{t)-l/a{t)^ J^ 



{f{t + r,t,x) - f{t,t,x)) m{dx) - g{t) 



0. 



(C14) 3Ku > such that, \/v eU,\fue U, 
1 



\v — u\ 



\f{v,v,x) - f{v,u,x)\ m{dx)<Ku. 



R 
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